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Abstract 

We consider the effect of CP violating phases in the MSSM on the relic density of the lightest 
supersymmetric particle (LSP). In particular, we hnd that the upper limits on the LSP mass 
are relaxed when phases in the MSSM are allowed to take non-zero values when the LSP 
is predominantly a gaugino (bino). Previous limits of mg < 250 GeV for < 0.25 can 
be relaxed to mg < 650 GeV. We also consider the additional constraints imposed by the 
neutron and electron electric dipole moments induced by these phases. Though there is 
some restriction on the phases, the bino mass may still be as large as ~ 350 GeV and certain 
phases can be arbitrarily large. 


It is well known that by considering the cosmological relic density of stable particles, 
one can establish mass limits on these particles. For example, light neutrinos {rrii, < 1 
MeV) contribute too much to the overall cosmological mass density if rrii, > 25 eV for 
< 1/4 [^. The relic density of heavier neutrinos are determined by the their annihilation 
cross section and yield lower bounds to neutrino masses rUi, > 4-7 GeV [@, §]• In the 
minimal supersymmetric standard model (MSSM), because of the many unknown parameters 
it becomes a significantly more complicated task to set limits on the mass of the lightest 
particle, the LSP, e.g. the annihilation cross section will depend on several parameters which 
determine the identity of the LSP. Here we will consider only neutralinos as the LSP (for a 
recent discussion on the possibilities for sneutrino dark matter, see @]). The mass limits will 
then depend sensitively on the masses of the scalars [^, and in general, limits are obtained 
on portions of the parameter space Q . For example when parameters are chosen so that the 
LSP is a gaugino (when the supersymmetry breaking gaugino masses Mi ^2 are taken to be 
smaller than the supersymmetric Higgs mixing mass, /i), for a given set of scalar masses, the 
requirement that Qh‘^ <1/4 places a lower bound on the gaugino mass. For small gaugino 
masses, this corresponds to a lower limit on the photino mass 0. At larger gaugino masses 
(but still smaller than p), the LSP is a bino 0. As the bino mass is increased, there is an 
upper limit mg < 250 GeV Hi. In p, it was shown that this upper limit is sensitive to 
the level of sfermion mixing. Here, we hnd that the upper limit to the bino mass is relaxed 
when phases in the MSSM, and in particular the phases associated with the off-diagonal 
sfermion masses, are allowed to take non-zero values. Now, the upper limit is increased to 
mg < 650 GeV. 

There has been a considerable amount of work concerning phases in the MSSM. For the 
most part these phases are ignored because they tend to induce large electric dipole moments 
for the neutron 0,0. To suppress the electric dipole moments, either large scalar masses 
(approaching 1 TeV) or small angles (of order 10“^, when all SUSY masses are of order 100 
GeV) are required. For the most part, the community has opted for the latter, though the 
possibilities for large phases was recently considered in |^. To reconcile large phases with 


small electric dipole moments, some of the sparticle masses are required to be heavy. In [13 


either large sfermion or neutralino masses (or both) were required. However, unless i?-parity 
is broken and the LSP is not stable, one would require that the sfermions be heavier than 
the neutralinos. If they are much heavier, this would result in an excessive relic density of 
neutralinos. The object of this letter is to determine the relationship between potentially 
large phases in the MSSM and the relic density while remaining consistent with experimental 
bounds on the electric dipole moments. 
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In the MSSM, the possibility for new phases arises from a number of sources. First, in 
the superpotential, there is the Higgs mixing mass, p. There are also several parameters 
associated with supersymmetry breaking: gaugino masses Mi, i = 1-3 (we will assume GUT 
conditions on these masses so that we need only consider one of them, M 2 ); and in the 
scalar sector, soft scalar masses; soft bilinear, Hp, and trilinear, Af, terms. Not all of these 
phases are physical |]T^. It is common to rotate away the phase of the gaugino masses, and 
to make Hp real, which ensures that the vacuum expectation values of the Higgs fields are 
real. We will also, for simplicity, ignore generation mixing in the sfermion sector though we 
will include left-right mixing. If furthermore we assume that all of the A/s are equal, and 
simply label them by A, we are left with two independent phases, the phase of A, 6a, and 
the phase of fi, 9^. The phase of B is hxed by the condition that Bfi is real. 

As both the relic cosmological density of neutralinos and the electric dipole moments are 
strongly dependent on the sfermion masses, it will be very useful to identify the combination 
of phases that enter the sfermion mass^ matrix. We take the general form of the sfermion 
mass^ matrix to be 


( 1 ) 


/ Ml + m'j + cos2j3(T3f — Qfsin^ 9w)Mz mfrUfA'^f \ 

\ Ml^ + m'j + cos 2f3Qf sirA OwM^ ) 

where are the soft supersymmetry breaking sfermion mass which we have assumed 

are generation independent and generation diagonal and hence real. We will also assume 
Ml = Mr. ruf is the mass of the fermion /, tan/3 is the ratio of Higgs vevs, and Rf = 
cot/3 (tan/3) for weak isospin -1-1/2 (-1/2) fermions. Due to our choice of phases, there is a 
non-trivial phase associated with the off-diagonal entries, which we denote by 
of the sfermion mass^ matrix, and 


m 


= Rffx + M = Rf\f^\R^^ + |A|e 




( 2 ) 


For a given value of /i and A, there are then two phases which can be distinguished (by Rf), 
and we denote them by 7 * and 'fh- The sfermion mass^ matrix is diagonalized by the unitary 
matrix U, 

/ u v\ 

. ( 3 ) 

\-v u 


where + \v\^ = 1, and we have taken u real. 

Previously we considered the effect of sfermion mixing on the relic density when the 
neutralino is mostly gaugino, and in particular a bino. The LSP is a bino whenever M 2 < fi 
and M 2 > 200 GeV; for smaller M 2 , the bino is still the LSP for large enough values of 
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yU. The bino portion of the M 2 — /i LSP parameter plane is attractive, as it offers the 
largest possibility for a signihcant relic density p!5| . The complementary portion of the 
parameter plane, with p < M 2 , only gives a sizable density in a limited region, due to the 
large annihilation cross sections to W^W~ and ZZ and due to co-annihilations with the 
next lightest neutralino (also a Higgsino in this case), which is nearly degenerate with the 
LSP We will also focus on binos as the LSP here. In addition to resulting in a sizable 
relic density, the analysis is simplihed by the fact that in the nearly pure bino region, the 
composition and mass of the LSP is not very sensitive to the new phases. However, as we 
will now show, the relic density, which is determined primarily by annihilations mediated by 
sfermion exchange, is quite sensitive to the phases, 7 * and 7 ?,. 

We begin by exploring the effect of the new phases on the relic density of binos. In 
the absence of these phases and in the absence of sfermion mixing, there is an upper limit 
0, ^ on the bino mass of mg < 250 GeV for VLh? <1/4. (This upper limit is somewhat 
dependent on the value of of the top quark mass. In |^, |^, we found that the upper limit 
was ~ 250 GeV for m^ ~ 100 GeV. When rrit = 174 GeV, the upper limit is 260 GeV. For 
rrit ~ 200 GeV, this limit is increased to mg < 300 GeV. Furthermore, one should be aware 
that there is an upward correction of about 15% when three-body hnal states are included 
n| which raises the bino mass limit to about 350 GeV. This latter correction would apply to 


the limits discussed below though it has not been included.) As the bino mass is increased, 
the sfermion masses, which must be larger than mg, are also increased, resulting in a smaller 
annihilation cross section and thus a higher relic density. At mg ~ 250 GeV, even when the 
sfermion masses are equal to the bino mass, ~ 1/4. Note that this provides an upper 
bound on the sfermion masses as well, since the mass of the lightest sfermion is equal to the 
mass of the bino, when the bino mass takes its maximum value. When sfermion mixing is 
included [Q, the limits, which now depend on the magnitude of the off-diagonal elements 
rufruf, are modihed. We hnd that this upper limit is relaxed considerably when the phases 
are allowed to take non-zero values. 

The dominant contribution to bino annihilation is due to sfermion exchange and is derived 
from the bino-fermion-sfermion interaction Lagrangian, 


= -j^g'iYRfPLBfR + YjPRBh)+h.c. 

= 75 fi '7 {fixPi + hwPi + fiyPR - hzPR) B + h.c. (4) 

where x = -Yrv*, y = Yru, w = Yru, z = -Ylv, Pr^l = (1 ± 75)/2, Yr = 2Qf and 
Yl = 2(Q/ — Tsf) where Qf is the fermion charge and T^f is the fermion weak isospin. 
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We compute the relic density by using the method described in ref. [||. We expand (avrei) 
in a Taylor expansion in powers of T/ 

{avrei) =a + b{T/ mg) + O (^(t/ mg)(5) 
The coefficients a and b are given by 


a — ^2 W®/ 
/ 


/ 


bf+1-3 + 


Sm'j 


( 6 ) 


(7) 


where df and bf are computed from the expansion of the matrix element squared in powers 
of p, the incoming bino momentum, and Vf = {1 — mj/m-)^/^ is a factor from the phase 
space integrals. 

We summarize the result by quoting the computed expression for df: 




ttf 


1287r 


Ai(m/ic^ — 2m^wz + mfz2 + /S. 2 {mfX^ + 2m-^xy + mfy^ 


A 1 A 2 


( 8 ) 


where Aj = m^ + m- — mj. The result for bf is too lengthy for presentation here, but was 
computed and used in the numerical integrations to obtain the relic densities. The results 
reduce to the results quoted in in the limit of zero phases. 

We show our results in Figures 1 and 2 for the upper and lower limits on mg as a 
function of the magnitude of the off-diagonal term in the top-squark mass^ matrix, m^, given 
the conditions 1) flgh^ <1/4, 2) the lightest sfermion is heavier than the bino, and 3) the 
lightest sfermion is heavier than 74GeV|]T^. In both hgures we have taken tan/? = 2 and 
mtop = 174 GeV. In Figure 1, |/r| = 3000 GeV and in Figure 2, |/i| = 1000 GeV. The various 
curves are labeled by the value of 7 ;, assumed, and in addition, mg has been maximized for 
all allowed values of 6^. The lower limit on m~ assumes 7 ;, = 7 r/ 2 . As one can see, when 7 ^ 
is allowed to take its maximal value of 7 r/ 2 , the upper limits are greatly relaxed to mg < 650 
GeV. With |p| and 7 b hxed, for a given value of m^ and 9^ all of the remaining quantities such 
as 1 ^ 1 , 0 ^, 7 b, and nif, are determined (though not necessarily uniquely, as some are double 
valued). 

The reason for the change in bounds becomes evident if we consider Ai A 2 = A and 


rUf << irt'. 


m 


5). Then 


«/ 


/\4 


(^0 

87rA2 


m- YI Im(u)^ -|- 0{mfm~) 
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( 9 ) 


) 2 A/^2 a^2 2 I i 2 * 2 i \ 

= I'^l + 0{mfm^) 

The size of u\v\ depends on the quantity r = mjlTif/{cos2(3{T^ — 2Q/sin^6*iy)M|); for 
r Ri 1, M |u| Ri .45, while for r -C 1 ,m |n| ~ r. Note that a non-zero value for siny^. removes 
the p-wave suppression for the fermion / (to a much greater extent than sfermion mixing 
alone |^) and greatly enhances the annihilation cross-section. Of course, for annihilations 
through the top quark channel, the p-wave suppression is not terribly strong, as mt is not 
much smaller than m-x. 

Jd 

For large m*, the diagonal mass terms must be taken large to ensure that the mass of 
the lightest stop is > mg. This drives up the masses of the other sfermions and suppresses 
their contribution to the annihilation. As rUt is decreased, must drop, the other sfermions 
begin to contribute and the upper bound on mg is increased. In particular, annihilation to 
p’s and r’s becomes important, since 


^ R 


ll,T 


^ R 


c,t 


W2w2 
■'L ^ R 


s,b 


= 81 


( 10 ) 


Decreasing reduces the effect of p’s and r’s, and this can be seen as a decrease in the upper 
bounds in Figures 1 and 2. For m^ sufficiently small, the stops become unmixed, diminishing 
somewhat their contribution and slightly decreasing the upper bound on mg. 

Lastly, we mention how the top contribution depends on y^. For the large values of mg 
allowed for rUt < 1500 GeV, top annihilation feels an enhancement from a non-zero sin y^ as in 
(§) (but note that for the top, the higher order terms in m^ are not small). However, except 
in a narrow region near rUt = 0, the two stop eigenstates ti and t 2 have a large mass splitting, 
and the bino couples with different strengths to ti and t 2 - Taking y^ away from tt decreases the 
coupling to the lighter stop (which can offer a much greater contribution to the annihilation 
than can the heavier stop), and we hnd that larger annihilation rates into tops are generally 
found nearer to y* = vr than y^ = 7r/2. In Figure 2, the sudden drop in the 9^ = n/S curve 
at 1500 GeV occurs because for m* > (tan (3 — cot /?) |p| , y^ — 7r/2 < y^ < y^ -|- 7r/2, and so 
for = 7r/8, yt is prevented from approaching vr. 

We turn now to the calculation of the electric dipole moments of the neutron and the 
electron. The EDM’s of the electron and quarks receive contributions from one-loop di¬ 
agrams involving the exchange of sfermions and either neutralinos, charginos, or (for the 
quarks) gluinos. In the case of the neutron EDM, there are additional operators besides the 
quark electric dipole operator, = \q(TiiuqF^'^ in which contribute. They are the gluonic 
operator Og = —\f°'^^GaGhGc HH and the quark color dipole operator, Oq = \qo'^^qT°‘G{(^ 
0. The gluonic operator is the smallest pT|, ^ when all mass scales are taken to be equal. 
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These three operators are conveniently compared to one another in and relative to the 
gluino exchange contribution to the operator, it is found that 0-y : Og : Og = 21 : 4.5 : 1. 

Because of the reduced importance of the additional operators contributing to the neutron 
EDM, we will only include the three contributions to the quark electric dipole moment. 
The necessary CP violation in these contributions comes from either 7 ^ in the sfermion 
mass matrices or in the neutralino and chargino mass matrices. Full expressions for the 
chargino, neutralino and gluino exchange contributions are found in []^|. The dependencies 
of the various contributions on the CP violating phases can be neatly summarized: for the 
chargino contribution 


dj ~ sin Q 


fj. 5 


with essentially no dependence on 7 /; whilst for the gluino contribution 

Trif smjf , 




( 11 ) 


( 12 ) 


independent of and the neutralino contribution has pieces that depend on both sin 6^ and 
WifSin'-yf. All three contributions can be important (including the neutralino contribution, 
in the case of the electron EDM), and depending on sin 6 *^ and sin 7 /, they can come in with 
either the same or opposite signs. In particular, sign[(ij/dj] = sign[sind^/siny^]. For the 
mass ranges we consider, the dipole moments fall as the sfermion masses are increased, and 
sfermion masses in the TeV range can bring these contributions to the neutron and electron 
electric dipole moments below the experimental bounds of \dn\ < 1-1 x 10 “^^e cm 


23 and 


|de| < 1.9 X cm [^, even for large values of the CP violating phases |^. However, 

these large sfermion masses are inconsistent with the cosmological bounds mentioned above, 
where sfermion masses must be relatively close to the bino mass in order to keep the relic 
density in check. 

We proceed as follows. We fix the value of 7 ;, and take \fi\ = 3000 GeV. Then for several 
values of fUt between 0 and 1500 GeV, we determine the upper bound on mg, as a function 
of 6fj,. As we vary across its full range, m^ and 7 * change, and this affects the annihilation 
rate (see (^ and following discussion), and consequently the bound on mg. Taking mg at 
its maximum value allows us to take as large as possible; although the electric dipole 
moments depend on mg as well, we hnd that the dependence on is sufficiently strong 
that the EDM’s take their minimum values for the maximum values of mg and M|. We then 
compute the quark and electron EDM’s as a function of and m^, and use the nonrelativistic 
quark model to relate the neutron EDM to the up and down-quark EDM’s via 


dn 2 (ddj^ du) . 


(13) 
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If we find no region of the 9^-Wit parameter space which satisfies both the neutron and 
electron EDM bounds, we decrease 7 f, and repeat the procedure. In practice, we hnd the 
bound on the neutron EDM the more difficult of the two to satisfy, and every region of 
the parameter space we show which produces an acceptable neutron EDM also produces a 
sufficiently small electron EDM. We will consequently drop further discussion of the electron 
EDM and concentrate on the neutron. 

For the large value of |/r| = 3000 GeV, we hnd that the largest contribution to the neutron 
EDM comes either from gluino exchange (for the more negative values of 9^) or chargino 
exchange (for the more positive values of 6*^), and that the value for \dn\ is too large unless y?, 
takes a relatively small value. In particular, we hnd non-negligible experimentally acceptable 
regions of the parameter space only for < 7r/25. In Figure (3) we show a contour plot of 
the neutron EDM as a function of 9^ and for = 7r/40. The shaded regions demarcate 
the range of for this choice of y?,. Much of this range produces a sufficiently small \dn\- 
As we increase mt, the d and u masses become large and \dn\ falls. As we move to values 
of TUt greater than ~ 1500 GeV, we begin to require a signihcant tuning of to produce 
< 1/4. 

Near the boundaries of the allowed range of 9^, y^ approaches jb ± 7r/2. As we explain 
above, the top contribution to the bino annihilation drops oh as we move away from y* = 
TT, and the upper bound on the bino mass, and sfermion masses, falls. It is these lower 
values for the sfermion masses which are primarily responsible for the sharp rise in \dn\ near 
the boundaries. There is one last subtlety which requires mention. For every set of the 
parameters {|/r|, rUt, y?,, 6*^}, there are two possible values for {yt,rn(,}. We hnd that one of 
the two sets of values always gives a smaller \dn\, given the cosmological constraints on the 
bino and sfermion masses, and it is these smaller values of \dn\ which we plot. 

We repeat this procedure for |/i| = 1000 GeV. For lower \fj,\, the chargino exchange 
contribution is enhanced relative to the gluino exchange contribution. In Figure (4), we 
show a contour plot of dn for y^ = n/8. At {fUt = 1500 GeV, 9^ = 0),dn vanishes, as 

= 9^ = sinyt = 0. Of course the gluonic and quark color dipole contributions to dn will 
not vanish everywhere along the contour dn = 0, but their contributions are < 3 x 10“^®ecm 
in the region plotted. Also, the gluino and chargino contributions scale diherently as mg and 
nij are changed. If we take mg and mj less than their maximum values in Figure (2), the 
contours of Figure (4) will shift, by < 1. x 10“^^. One should therefore concentrate on the 
qualitative features of Figure (4), as the exact positions of the contours are not signihcant. 

In summary, we have found that CP violating phases in the MSSM can signihcantly 
affect the cosmological upper bound on the mass of an LSP bino. In particular, taking the 
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maximal value 7r/2 for the phase 'jb of the off-diagonal component of the T 3 = —1/2 sfermion 
mass matrices pushes the upper bound on mg up from ~ 250 GeV to ~ 650 GeV. When we 
additionally consider constraints on neutron and electron electric dipole moments, we hnd 
the upper bound on mg is reduced to ~ 350 GeV. Various combinations of the CP violating 
phases are constrained as well; \9^\ < 0.3 and | 7 fe| < n/6 for |/r| > 1000 GeV, while ■jt and 
6a are essentially unconstrained. We note that although the bounds on and gf, are small, 
they are much larger than the values of order 10 “^ typically considered. 
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Figure Caption 


Fig. 1) Upper limits on the bino mass as a function of the off- diagonal element frit 
in the top squark mass^ matrix, for various values of 'jb, the argument of the 
off-diagonal element of the T 3 = — 1/2 sfermion mass^ matrix. Also shown is 
the lower bound (lowest curve) on the bino mass assuming 7 ^ = 7 r/ 2 . The 
value of |/i| was chosen to be 3000 GeV. 

Fig. 2) As in Fig. 1, with \n\ = 1000 GeV. 

Fig. 3) Gontours of the neutron electric dipole moment, dn, in the - frit plane in 
units of 10“^®e cm. The value of |/i| was chosen to be 3000 GeV. The shaded 
region corresponds to values of 6 *^ and frit which are not allowed algebraically 
for this value of /i and 7 ;,. 

Fig. 4) As in Fig. 3, with \fi\ = 1000 GeV. 
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